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We aim to offer a kind of unifying view on two popular topics in the studies of nonperturbative aspects of 
Yang-Mills theories in the Landau gauge: the so-called Gribov-Zwanziger approach and the Kugo-Ojima con- 
finement criterion. Borrowing results from statistical thermodynamics, we show that imposing the Kugo-Ojima 
confinement criterion as a boundary condition leads to a modified yet renormalizable partition function. We 
verify that the resulting partition function is equivalent with the one obtained by Gribov and Zwanziger, which 
restricts the domain of integration in the path integral within the first Gribov horizon. The construction of an 
action implementing a boundary condition allows one to discuss the symmetries of the system in the presence 
of the boundary. In particular, the conventional BRST symmetry is softly broken. 



I. INTRODUCTION 

The Gribov-Zwanziger (GZ) approach focuses on the is- 
sue of gauge copies in the Landau gauge. Gribov signalled 
in his seminal work |1] that the Landau gauge condition, 
d^Afj = is ambiguous: there exist gauge equivalent config- 
urations A'^ which also obey 3^A^ — 0. Examples of gauge 
copies are provided by the zero modes of the Faddeev-Popov 
(FP) operator, which enters the quantization formula of Yang- 
Mills theories. Indeed, given an infinitesimal gauge transfor- 
mation connecting with A^, i.e. A^ — A", — D"l'(i)'\ it is 
clear that 8^A^ = d^A^ = is fulfilled when M"*©* = 0, with 

M"'' = -d^,Df = -d^ (3^8"'' + 8f"''A'jj) being the FP oper- 
ator. We recall that the FP action in the Landau gauge for a 
J-dimensional Euclidean gauge theory, with d <4, reads 

5YM+gf = Sym + I d'x (b"d,A; + l^'d.Dfc'^ , ( 1 ) 

with 5ym = \ j d^'xFl^yjF^yj the classical Yang-Mills action. 
Expression ^ enjoys the well-known BRST symmetry, gen- 
erated by the nilpotent operator s, ~ 0, i.e. 

sAl = -Dfc'' , sc" = ^gf"'"'c''c' , sc" = b" , sb" = . (2) 

For the partition function, we can write 

Z = Jd^tFP = J d^e-^^^+i^' = j dAA&iM"''h{dA)e-^^^ . 

(3) 

We introduced the notational shorthand <I> denoting all the 
fields present in the action, with dp-pv the usual FP measure. 

Gribov proposed to restrict the domain of integration to 
the subspace Q., where the Hermitian operator M"^ is posi- 
tive definite. More precisely, we define the Gribov region as 



D. = {A"i, 3^,A", = 0, M"* > 0}. We recognize that configura- 
tions A'^EQ.aie relative minima of the functional / d'^x{A'^)^, 
u e SU{N). The boundary, dO., of il is called the (first) Gribov 
horizon. It was shown with increasing rigor that Q. is convex, 
bounded in all directions in field space, and that every gauge 
field has at least one gauge equivalent representant in fi (see 
12, [H and references therein). The inverse of the FP operator, 
or equivalently the ghost propagator with external gauge field, 
G"^{k^A), can be used to implement the restriction to fi, as 
done semiclassically by Gribov. Following [1], we can write 

At lowest order, it can be shown that 1 + 0(^, A) is a decreas- 
ing function of ^ Jll], hence one can impose 

l+a(0,A)>0. (5) 

Condition (|5]l, known as the Gribov no-pole condition, implies 
that the ghost propagator G"^{k,A) has no poles at finite non- 
vanishing k. Moreover, positivity of G"^{k,A) ensures that the 
Gribov horizon dQ. is not crossed. As done by Gribov Ol], the 
no-pole condition can be embodied into the partition function 
using a 5-function\ 

Z' = y^c/4>8(l+a(0,A))e"'^™+8f. (6) 

Later on, Zwanziger 120 was able to implement the no-pole 
condition to all orders. Relying on the equivalence between 
the microcanonical and the canonical Boltzmann ensemble 
(see also Section II), he was able to show that the partition 
function (|6]l has to be replaced by 

Z" ^ / t/A5(a^,A^,)detM"''e-^™+//rf''.^'>W ^ (7) 



' Condition ^5) can be implemented by inserting a step function factor 
0{l + <s(k,A)). However, in the thermodynamic limit, the 0-function can 
'Electronic address: |david.dudal@ugent.be, sorella@uerj.br, nele.vandersickel@ugeilld)e^tBri'fli5irajditli*p®fi^oiilctiDii, see rn.l3. rt5n . 
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where the Zwanziger horizon function reads 

I d'xh{x,A)^g^l d''xd''yf'%{x){M-'r%x,y)f^'Al{y) 

(8) 

The mass parameter y is determined by a gap equation, com- 
monly called the horizon condition 

{h{x))^d[N^~\), (9) 

where {h) is calculated with Z", i.e. with the measure 
d/jFpe'^!''''''^^''\ The factor £/(A?2 _ 1) jjjtljgj-lj g was ob- 
tained by determining the lowest eigenvalue of the FP op- 
erator Working out the condition (|9]l at lowest order repro- 
duces the Gribov result 1 1 ] . The action corresponding to the 
partition function O contains the nonlocal horizon term 
To arrive at a workable quantum model, it was shown that 
^ can be put in an equivalent local form by introducing a 
set of complex conjugate commuting variables, (9^' , cp^',' ), and 
anticommuting ones, (co",', co^' ), so that we finally obtain the 
Gribov-Zwanziger action, 

sgz = 5YM+gf+/rf''x((r;'3vZ)fc-co;;'3vZ)f< 

-g{dvTo";)f'"^D^c)\;^ 

-I'gf'K (<■+<■) +d{N^ - 1)y4) . (10) 

The horizon condition (|9]l is translated as ^ = 0, with r(Y) 

the effective action, defined as e^^ — J d'Pe^^'^^. This can be 
easily checked, given that we take y 7^ 0. The mass parameter 
Y turns out to be proportional to A^s, and as such it can give 
rise to nonperturbative corrections. This is not unexpected, 
as the restriction to the region i2 is a highly nontrivial oper- 
ation, which goes beyond perturbation theory. At the pertur- 
bative level, the ghost propagator stays positive. We are thus 
far from the horizon and nothing happens. It is only at lower 
momenta, where normal perturbation theory starts to fail, that 
the fields begin to feel the restriction to Q.. Having brought the 
action in standard local form, we have all the usual concepts 
and machinery of local quantum field theory to our disposal. 
A first important property of ( fTOb is its renormalizability to 
all orders of perturbation theory. Hence, the restriction to Q. 
makes perfect sense at the quantum level, and finite results are 
found, consistent with the renormalization group |3, 4]. We 
stress here that the action ( fTOl ). with the horizon condition (|9]l 
implemented, is nothing else than the correct extension to all 
orders of the usual Yang-Mills action, supplemented with the 
Landau gauge fixing, in the presence of a nontrivial boundary 
condition, being the no-pole condition (|5]). In this fashion, it is 
assured that we have taken care of a certain amount of gauge 
copies, including those related to the zero modes of M"^. No- 
tice that this does not mean that the Gribov issue has been 
completely solved. It is known that £2 still contains copies, 
related to the fact that / d''x{A'p^ can have many relative min- 
ima starting from the same A^. A further restriction is needed, 
keeping only gauge configurations that are absolute minima of 
/ d''x{A^)^; the latter define the fundamental modular region 
(FMR) A. 



Evidently, the extra fields, (9^'', cp^jj'', can influence 

the dynamics of the theory in a nontrivial fashion [5]. These 
fields arise as a consequence of the presence of the Gribov 
horizon. As such, they can give rise to additional nonperturba- 
tive effects. For example, in L^, we have provided evidence of 
the existence of a dimension 2 condensate, ((p^'^cp^'' - ^'^CO^'^), 
md = 4. A posteriori, this is not that surprising, given that the 
restriction to Q. introduces the mass scale yinto the theory, and 
that the horizon condition (|9]) can be reexpressed at the local 
level as (g/"^'A;',((p;[;'' + (p^')) = -2fd{N^ - 1), i.e. a dimen- 
sion 2 condensate for d = 4. Nontrivial condensates are an 
important source of nonperturbative effects in gauge theories, 
hence the general interest in their study. In particular, dimen- 
sion 2 condensates attracted a lot of attention in recent years, 
see e.g. and references therein. In the current case, the 

operator (9^' 9^' — co^'^^co"' ) can be added to the theory in a 
way that preserves renormalizability J^, which is already a 
remarkable feature, indicative of its possible relevance. We 
studied the effects of this condensate using variational pertur- 
bation theory, and found that the gluon propagator^ does not 
vanish at zero momentum (D{Q) ^ 0), that the ghost propa- 
gator behaves like ~ p- at small momenta, and that there is 
a violation of positivity in the gluon propagator |6]. Any of 
these findings is in good agreement with all most recent lattice 
data, obtained at previously unseen large volumes [8, 9]. Also 
certain results based on Schwinger-Dyson (SD) and/or Func- 
tional Renormalization Group (FRG) equations are consistent 
with these data, see e.g. ifTol [Till . Without taking into ac- 
count the effects related to '^"^^^^ -'^'^^''^), the GZ action 
( [Tol l also leads to the positivity violation of the gluon propaga- 
tor, however with Z)(0) = 0, and an infrared enhanced ghost. 
These latter two results are no longer supported by lattice data. 
Hence, it seems crucial to take into account additional nonper- 
turbative effects related to the restriction to the region £2 (i.e. 
the boundary condition) to allow for consistency between the 
analytical GZ results and most recent lattice predictions. The 
interpretation of the analytical |5, 6] and the lattice results of 
[|8|] was challenged in papers like [9, 10, J2J. It was argued 
that the ghost propagator must be infrared enhanced to en- 
sure confinement, whereby only colorless states are physical. 
These statements are based on the Kugo-Ojima (KO) analysis 
of gauge theories lfl3[fl4ll . This analysis relies on the operator 
formalism, and it has been shown that, given a globally well- 
defined BRST charge Qb, the color charge is a BRST exact 
variation, = Qb{. . .), if the gluon propagator contains no 
massless poles. The color charge Q" is then well-defined only 
if the KO confinement criterion holds 

m(0) = -1, (11) 

with u{k^) defined through the following Green function 

I Ae'*^^ (^D";'c''ix)D'"''c%0))^^ = 5"^ (^Pi,v{k)u{k^) - ^ 

(12) 



- The Landau gluon propagator can be parametrized in terms of the form 
factor D{k^) as {Af,{k)A'i{-k)) = S''''{8,„ - '^)D(k^). 
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in Minkowski space. (. . .)pp stands for the expectation value 
taken with the FP action ([5]l, while P/jvik) = g/jv ^ k^k^/k^ 
for the transverse projector. Using the nilpotent BRST charge 
Qb, one can invoke its cohomology to define the physical sub- 
space, and by means of Q" = Qb{---), conclude that phys- 
ical states cannot carry color A few comments are in or- 
der First of all, in the KO framework T4], the existence 
of a globally well-defined BRST charge is assumed. Thus, 
the issue of the (non)existence of a nonperturbatively valid 
BRST symmetry is not explicitly faced. Secondly, the link 
between the BRST charge Qg and global color charge Q" is 
made using the action ([T]), i.e. by employing the usual FP 
gauge fixed action. As such, the Gribov problem is simply not 
addressed. It is worth noticing that Kugo and Ojima did not 
impose the criterion ( fTTT i. but they derived it as a condition 
to be checked/calculated. Though, nowadays, in functional 
formalisms as in ifTol] . the criterion is used as input. Kugo 
showed in ifTsll that, in the Landau gauge, one can rewrite the 
ghost propagator 



^"''^^^ k^ \+u{k^)+k'^v{k^) ' ^^^^ 

meaning that the criterion (fTTT l is equivalent to an infrared en- 
hanced ghost. The ghost enhancement is then imposed as 
a boundary condition in order to favor the so-called scaling 
type solution of the SD and/or FRG equations |10]. Let us 
already draw attention to the close similarity existing between 
the no-pole condition (|5]l and the criterion ( fTTT l. Imposing^ 
a(0) = — 1 exactly corresponds to m(0) = — 1. 



II. m(0) = 1 AS A BOUNDARY CONDITION 

We want to show that the constraint m(0) = —1 can be im- 
plemented directly into the theory, by appropriately modifying 
the measure one starts from. We shall see that the resulting ac- 
tion will be exactly the same as the GZ action. This has several 
interesting consequences which we will discuss in Section III. 
We shall first give an overview of some results from thermo- 
dynamics we intend to employ. 



A. Microcanonical ensemble and equivalence with the 
canonical Boltzmann ensemble in the thermodynamic limit 

We consider a discrete system, whose Hamiltonian is 
H{q,p), with 3A^ degrees of freedom. The averages in the 
microcanonical ensemble are constructed out of 

E(£)= / dn = ( dnb{E-H), 
Jh=e J 

where d/j — d^^qd^^p represents the classical phase space 
and E stands for the constant energy of the system. Aver- 
ages in the microcanonical ensemble are defined by (0)^^^^.^ = 



^ 0(0) is related to c(0,A) by making the gauge field dynamical and per- 
forming the coiTesponding path integration. 



^^Y^^-j^- In order to establish the equivalence between the 
microcanonical and the (Boltzmann) canonical ensemble we 
rewrite the quantity in the following form 

E(£) = fdf,8{E-H)= f d^i ^e^iE-H) 
J J J-ioo+t 2ni 



/(P)= I d^ie 



(14) 

co(p)--log/(p). (15) 



It can be shown that, in the thermodynamic limit, N,V ^ °°, 
with N/V fixed, the saddle point approximation becomes ex- 
act. We refer to Iil5il for an overview of the proof. So, 



E(^) = ^/(n,whhco'(n = ^-o. 



Jd^iHe-^*" 



From eq. (fTSI) it follows that 



(H) 



Boltz 



JdlJe 



(16) 



(17) 



This is the gap equation determining the critical parameter 
P*. Analogously, it can also be shown that IJ5.1 {0)f^^^^ = 

('5>Boitz = '^jX^,-p"«" for average of any quantity 0{q,p). 



B. Imposing the KO criterion yields the GZ framework 

Starting from (fT2l i and performing Lorentz and color con- 
tractions and taking the p ^0 limit, we can write 



-{VT)-' / d'y / d'x{D";'{x)D';;{y){M-'y^ix,y)) 



FP 



= iN^-l){{d-l)u{0)-l), 



(18) 



after passing to Euclidean space, as in any functional or lattice 
approach. VT denotes the spacetime volume. The identifica- 
tion between (. . . c'^' {x)c'' (y)) pp and (. . . {M^^Y'' {x,y))pp can 
be easily proven using the path integral (O. After discard- 
ing terms which are total derivatives, one easily sees that the 
quantity in the l.h.s. of (fTSI l is, up to the sign, the Zwanziger 
horizon function h{x). More precisely, we have 

^d''x{gf"'Al{x){M-'r{xM'""'A';{y))^^ 



-{VT)-' j d''x{h{x))pp^-{h)pp. 



(19) 



We observe that the KO condition cannot be realized with the 
standard FP measure li/jpF, otherwise we would have 



{h{x))pp = d{N^-\), 



(20) 



which would contradict Zwanziger's result (|9|l, obtained by 
restricting the path integral to the Gribov region £1. We now 
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implement the KO criterion m(0) = — 1 as a boundary condi- 
tion, amounting to start from the modified measure 

d/UFP d/j' = d/jpp 5 (vTd{N^ ^^^^ J ^''■^'^(■^)) ' "^^l) 

which clearly implements {h{x)) = d{N'^ — 1), or equivalently 
m(0) — —1. We are thus led to consider the partition function 



Jd/j' = J d/jFP 8 \VTd{N^ ^^^^ j '^''^^^^ 

dA ?>{dA) detfAf c^^ym 5 (vTd{N^ J ^''-"^^^ 



/' 



d<P5 



(vTd{N^ -1)- j d'^xh{x) 



,-5YM+gf 



(22) 



Expression ()22)) defines a microcanonical ensemble. Since we 
are working in a continuum field theory, we are working in 
the thermodynamic limit, hence we have an equivalence with 
a Boltzmann canonical ensemble as outlined in the previous 
section. Using analogous arguments as there, we arrive at 



where the mass parameter y follows from the gap equation 



d{N^-l) = {h{x)) 



Boltz 



/ d^lppe-'^^ h{x) 



(24) 



which is the analogue of ( fTTl ). We conclude that we can con- 
sistently encode the boundary condition ( fTTl ) at the level of 
the action, which turns out to be identical to the GZ action, 
eq.(|7]). Of course, we can localize it into the form (fTOl i. with 
corresponding local formulation of the gap equation. 



III. DISCUSSION 

Naively, one might already expect that the introduction of 
a nontrivial boundary condition can seriously influence the 
dynamics of the theory. One of our main points is to stress 
that one should introduce the boundary condition into the the- 
ory from the beginning, to fully grasp all its nontrivial as- 
pects. Having at our disposal an action automatically imple- 
menting the boundary condition, we can study an important 
aspect: the symmetries of the theory in the presence of the 
boundary. In principle, imposing a boundary could jeopar- 
dize certain symmetries of the original action. We have al- 
ready shown in |0] that placing a boundary in field space 
at the first Gribov horizon breaks the conventional BRST 
symmetry (|2|i. The practical implementation of the horizon 
by means of the GZ formulation confirms this, as sSgz — 
gf J d^'xf"'"^ (a;;©^'' - (Z)^™c'") (cp^;'' + (p^'^) ) O. We no- 
tice that the BRST generator (|2|i has a natural extension to the 
extra fields ((p^'',(p|[;'',co^'',a)|[;' ), given by 



forming 2 pairs of BRST doublets. In the absence of the GZ 
restriction or equivalently of the KO criterion, i.e. when y = 0, 
we are then assured that these fields are trivial in the BRST 
cohomology, thus completely decoupling from the physical 
subspace [16]. Let us come to another important statement. 
If Gribov copies are taken into account a la GZ, which is 
equivalent to imposing the KO criterion as we have verified in 
the previous section, the precise meaning of the KO confine- 
ment criterion becomes unclear. Since the BRST symmetry is 
broken, one can no longer simply use it to define the physi- 
cal subspace. It is sometimes mentioned in the literature that 
there might be a nonperturbative, globally well-defined BRST 
charge g^, and it is this Q'^ KO is referring to fMIUl. We 
cannot exclude this possibility, but this is a highly nontrivial 
statement and, obviously, it asks for a proof. At present, we 
are unaware of any such proof. Even if the charge 2b would 
be known, the KO analysis would need to be reworked from 
the start, as it explicitly relies on the FP action ([T]l and conven- 
tional BRST symmetry Simply stating that Q'g must exist 
in order to define the physical subspace analogously to what is 
done at perturbative level, does, in our opinion, not solve the 
problem. Also, the relation between a new BRST charge Qg 
and the global color charge would need to be reestablished, 
if any relation exists to begin with. One can speculate that it 
might be possible to modify s into Sy, such that limy-^o Sy = s 
and ^Y^Gz = 0. However, such a possibility can be easily dis- 
proved. Indeed, as y has mass dimension 1, and by keeping in 
mind that the BRST generator s does not affect the dimension 
of the fields'*, it is impossible to introduce extra y-dependent 
terms in the BRST transformation of the fields while preserv- 
ing locality, Lorentz covariance and global SU (N) structure. 
Let us briefly return to the functional SD (ERG) approaches. 
Now that a renormalizable action has been constructed, which 
implements the desired boundary condition explicitly, one can 
write down the corresponding SD (ERG) equations and try to 
solve them, given that the gap equation ( l24l l must be solved 
simultaneously. We expect that different kinds of solutions, 
similar to those found in [10 1, will emerge. A way to dis- 
tinguish between them could be based on selecting the most 
stable solution, i.e. the one with the lowest corresponding 
vacuum energy. We notice that there is still a lot of infor- 
mation available about the action ( fTOb . e.g. nonrenormaliza- 
tion properties typical of the Landau gauge, a renormalizable 
softly broken Slavnov-Taylor identity, etc. [6]. We conclude 
that in the current spirit of using the KO condition ( fTTl ) as in 
M\M, there is no clear connection between the KO criterion 
m(0) = — 1 and the highly nontrivial issue of confinement. All 
that one can say is that there is a violation of positivity in 
the gluon propagator, which is indicative of confinement, but 
certainly not a proof of it. Also, in the light of our previous re- 
sults, we disagree with the statement made in |l9l [Tol[T2|[ about 
the fact that the SD(ERG) solution with an infrared enhanced 
ghost propagator would refer to the absolute Landau gauge, 
i.e. to the restriction to the EMR A. Unfortunately, at present. 



(25) 



The usual canonical dimensions are assigned to the fields [It 
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a way to implement the restriction to the FMR A remains com- 
pletely unknown. Moreover, we remind that the recent lattice 
data have given quite clear evidence about the fact that the 
ghost propagator is not enhanced in the infrared, within the 
current accuracy of implementing the Landau gauge as the 
mmimum of the functional /t/''x(Ap^ u e SU{N) [8]. Even 
if in the future more powerful algorithms would bring the sim- 
ulations closer to the FMR A, the ghost propagator will not get 
more enhanced than before, on the contrary ilTll . Moreover, 
we have shown in this letter that the KO boundary condition 
is equivalent with the GZ framework, which explicitly refers 
to the restriction to the Gribov region Q.. This is irrespective 
of the fact that the ghost is enhanced or not, implementing KO 
breaks the conventional BRST symmetry, and refers to £2, not 
to A. We wish to underline that implementing the boundary 
as in eqs.(l23]l and ( l24b will not necessarily give rise to an in- 
frared enhanced ghost. Additional nontrivial quantum effects 
can combine with the boundary effect, we refer for instance to 
the effects of the operator cp^',' cp^' — W^w/^f in the case under 
study. We can make the analogy with spontaneous symme- 
try breaking: although the starting action in that case enjoys a 
certain symmetry, nonperturbative quantum effects can induce 
a shift from the symmetric, but unstable, vacuum^, causing 
qualitative and quantitative effects in the theory. 

In conclusion, we hope that this letter has clarified the re- 
lation between the KO and GZ framework. Our main result 
is expressed by eqs. (l23] l and (l24l l. which show that the KO 

^ We refer to condensates which are dynamically favored by lowering the 
vacuum energy. 



and GZ frameworks are equivalent, provided the KO bound- 
ary condition is properly taken into account from the begin- 
ning. The conventional BRST operator ^ suffers from a soft 
breaking, which relies precisely on the implementation of the 
boundary condition. Some ingredients in certain formalisms, 
which we have tried to outline, have thus to be considered 
as assumptions rather than as proofs. In particular, the pre- 
cise relation between implementing the KO criterion and con- 
finement remains to be clarified. One of the challenges lying 
ahead is how to define what the relevant physical operators are 
in the KOGZ framework, if there is no (local) nilpotent BRST 
symmetry generator found. 
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